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We outline the key-steps towards the construction of a physical, fully relativistic cosmology, in 
which we aim to trace Dark Energy and Dark Matter back to physical properties of space. The 
influence of inhomogeneities on the effective evolution history of the Universe is encoded in backre- 
action terms and expressed through spatially averaged geometrical invariants. These are absent and 
interpreted as missing dark fundamental sources in the standard model. In the inhomogeneous case 
they can be interpreted as energies of an emerging scalar field (the morphon). These averaged in- 
variants vanish for a homogeneous geometry, where the morphon is in an unstable equilibrium state. 
If this state is perturbed, the morphon can act as a classical inflaton in the Early Universe and its 
de-balanced energies can mimic the dark sources in the Late Universe, depending on spatial scale as 
Dark Energy or as Dark Matter, respectively. We lay down a line of arguments that is qualitatively 
conclusive, and we outline open problems of quantitative nature, related to the interpretation of 
observations. 
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I. GENERAL RELATIVITY AND COSMOLOGY 
A. The foliation issue and the notion of an effective cosmology 

The homogeneous-isotropic standard model of cosmology, being itself a particular solution of Einstein's general 
theory of relativity, does by far not exploit the degrees of freedom inherent in the geometry as a dynamical variable. 
It is this richer tone of general relativity - as compared to the Newtonian theory - that opens the possibility to 
generalize cosmological models, notably by including inhomogeneous structure also in the geometrical variables. There 
are several guidelines to be emphasized in such a generalization: firstly, a cosmology is thought of as an evolving space 
section that implies the need to speak of a foliated space time, introducing four degrees of freedom (the lapse and shift 
functions in an ADM setting). This necessarily implies, on general grounds, a breaking of four-dimensional covariance. 
This fact should not be confused with coordinate- or gauge-dependence of the resulting cosmological equations and 
variables, however. Secondly, a cosmology purports an effective point of view in the sense that the evolving spatially 
inhomogeneous variables are thought of as being "averaged over" in a way that has to be specified. We aim at a 
description that only implicitly refers to a metric. However, if a metric is to be specified, a cosmological metric is then 
to be considered as an effective, "smoothed out" or template metric, being not necessarily a solution of the equations 
of general relativity. Finally, a physical cosmology should be characterized by such an effective evolution model, 
an effective metric to provide the distance scale for the interpretation of observations, or alternatively an evolution 
model for average characteristics on the light cone, together with a set of initial data. These latter are to be related 
to physical properties of fundamental sources, but also to the geometrical data at some initial time (effective, i.e 
"averaged" quantities of known energy sources, intrinsic and extrinsic curvature). This latter clearly emphasizes the 
absence of any phenomenological parameters. Those would just parametrize our physical ignorance. All these points 
will be made explicit in what follows. 

B. The dark side of the standard model: postulated sources and proposed solutions 

The high level of idealization of the geometrical properties of space in the standard model leads to the need of 
postulating sources that would generate "on average" a strictly, i.e. globally and locally, homogeneous geometry. It 
is here where a considerable price has to be paid for a model geometry that obviously is not enough to meet physical 
reality, unless we really believe that we can find the missing sources: 96 percent of the energy content is missing 
in the form of a) a postulated source acting attractive like matter, so-called Dark Matter (= 23 percent) and b) a 
postulated source acting repulsive, so-called Dark Energy (= 73 percent). Evidence for the former does indeed come 
from various scales (galaxy halos, clusters and cosmological, see e.g. |63fl). while evidence for the latter only comes 
from the apparent magnitude of distant supernovae (see (30 . [ssl . [43 | for the latest data) that, if interpreted within 
standard model distances, would need an accelerating model. In the simplest case this volume acceration is achieved 
by a homogeneous-isotropic cosmology with a cosmological constant. It should be emphasized that when we speak of 
evidence, we already approach this evidence with model priors [36l. [68l. [70| . Keeping this idealization for the geometry 
of the cosmological model for example, one has to conjecture fundamental fields in proportion to the missing dark 
components on cosmological scales. The search for these fields is one major research direction in modern cosmology. 

Another huge effort is directed towards a generalization of the underlying theory of gravitation. While this would 
generalize the geometry of the model, it is not clear why all these efforts go into a generalization of general relativity 
and not into the generalization of the cosmological model within general relativity. There are certainly good lines of 
arguments and various motivations in particle physics and quantum gravity to go beyond the theory of Einstein (for 
reviews see [l^, [i^), but the "dark problem" may be first a classical one. 

Looking at generalizations of the standard model within general relativity can be identified as a third research 
direction to which we dedicate our attention here. In light of current efforts it is to be considered conservative, since 
it does not postulate new fundamental fields and it does not abandon a well-tested theory of gravitation , [I^ , 
[39| (for reviews on the physical basis of this third approach see d, [3 and [131 )• Among the works in this latter field, 
research that analyzes spherically symmetric exact solutions has been meanwhile developed to some depth, and has 
determined the constraints, necessary to explain Dark Energy, on a postulated observer's position within a large-scale 
void (see d, i, [H, ^ and references therein). 
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C. Fictitious and physical backgrounds 

Perhaps a reason for not questioning the standard model geometry within general relativity and to go for the 
search for fundamental fields or for generalizations of the laws of gravitation is the following belief: effectively, i.e. 
"on average", the model geometry has to be homogeneous, since structures should be "averaged over". Then, due to 
observational facts on large scales (the high degree of isotropy of the Cosmic Microwave Background, if the dipole is 
completely eliminated due to our proper motion with respect to an idealized exactly isotropic light sphere), and first 
principle priors (the strong cosmological principle that requires the universe model to look the same in all directions), 
the model geometry is taken to be locally isotropic. 

Taking this reasoning at face value we must note two points: the notions of homogeneity and isotropy in the 
standard model are both too strong to be realistic: firstly, local isotropy implies a model that is locally and globally 
homogeneous, i.e. despite the conjecture that the homogeneous model describes the inhomogeneous Universe "on 
average", this strict homogeneity does not account for the fact that any averaging procedure, in one way or another, 
would introduce a scale-dependence of the averaged (homogeneous) variables |28| . This scale-dependence, inherent in 
any physical averages, is suppressed. Even if a large scale of homogeneity exists (we may call this weak homogeneity 
principle), the model is in general scale-dependent on scales below this homogeneity scale j69| . The same is true for 
isotropy: while the averaged model may be highly isotropic on large scales, a realistic distribution on smaller scales 
is certainly not (we may call this weak isotropy principle). Correspondingly, a weak cosmological principle would be 
enough to cover the reality needs while still facing observational evidence on large scales. 

We may summarize the above thoughts by noting that, on large scales, a homogeneous-(almost)isotropic state 
does not necessarily correspond to a homogeneous-(almost)isotropic solution of Einstein's equations. These former 
states are the averages over fluctuating fields and it is only to be expected that the state coincides with a strictly 
homogeneous solution in the case of absence of fluctuations. In other words, looking at fluctuations first requires 
to establish the average distribution. Only then the notion of a background makes physical sense [4l|. Current 
cosmological structure formation models, perturbation theories or N-body simulations, are constructed such that the 
average vanishes on the background of a homogeneous-isotropic solution 01 • A such chosen reference background 
may be a fictitious background, since it arises by construction rather than derivation. On the contrary, a physical 
background is one that corresponds to the average (whose technical implementation has to be specified, and which 
is nontrivial if tensorial quantities like the geometry have to be "averaged"). A sound implementation of a physical 
background will be a statistical background where not only solutions but ensembles of solutions are averaged. Having 
specified such an averaging procedure, a physical cosmological model may then be defined as an evolution model 
for the average distribution. Despite these remarks it is of course possible that the homogeneous solution forms 
at the same time the average. A well-known example is Newtonian cosmology Q. It is also conceivable that the 
homogeneous solution provides, in some spatial and temporal regimes, a good approximation for the average. Still, it 
is important to consider perturbations on the correct background solution (40| . 

II. SCALAR FIELD MODELS AND THE MORPHON 

A. Effective evolution of inhomogeneous universe models 

Taking the point of view of generalizing the cosmological model within general relativity by abandoning the strong 
cosmological principle (strict homogeneity and isotropy on all scales) and replacing it by the weak cosmological 
principle (existence of a homogeneity scale and restriction to effective states that are almost isotropic on the scale 
of homogeneity) leads us to a "rewriting of the rules" to build the cosmological model. We shall consider the rules 
that led to the standard model of cosmology and replace them by their more general counterparts. It follows a 
basically similar framework that displays, however, a signature of inhomogeneity through the occurence of so-called 
backreaction terms and through a manifest scale-dependency. We shall not introduce new principles or assumptions, 
apart from the above outlined relaxation of the cosmological principle. We shall restrict ourselves to the simplest case 
of an irrotational dust model (for generalizations of the dust model Q with non-constant lapse function see [loj , and 
for additionally non- vanishing shift see [5., &, .3 3, , i44j ) . 

• As in the standard model we introduce a foliation of space time into fiow-orthogonal hypersurfaces. We generalize 
the notion of Fundamental Observers to those that are in free fall also in the general space time. Although, as 
in the standard model, this setup depends on the chosen foliation, we presume that this choice is unique as it 
prefers the fundamental observers against observers that may be accelerated with respect to the hypersurfaces. A 
general inhomogeneous hypersurface - contrary to the homogeneous case - will, in this setting, unavoidably run into 
singularities in the course of evolution. This is to be expected in a given range of spatial and temporal scales, since 
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we are treating the matter model as dust. This is not a problem of the chosen foliation, but a problem of the matter 
model that has to be generalized, if small-scale structure formation has to remain regular, and this can be achieved 
by the inclusion of velocity dispersion and vorticity. 

• As in the standard model we confine ourselves to scalar quantities. We replace, however, the homogeneous quantities 
by their spatial averages, e.g. the homogeneous density gnit) is replaced by {g)-p (t) for the inhomogeneous density 
g that is volume-averaged over some compact domain V. We realize the averaging operation by a mass-preserving 
Riemannian volume average. In some mathematical disciplines and in statistical averages at one instant of time, it 
may be more convenient to introduce a volume-preserving averager, but thinking of an averaging domain that is as 
large as the homogeneity scale we have to preserve mass rather than volume. Furthermore, the average is performed 
with respect to the above-defined Fundamental Observers. Spatially averaging a scalar ^/{t,X^), as a function of 
Gaussian coordinates X* and a synchronizing time t, is defined as: 

{^{t,X^))v:^^ ^{t,X^)d^ig , (1) 

with the Riemannian volume element d^g :— ^Jgd^X, g := det((7ij), and the volume of an arbitrary compact domain, 
VT,{t) := /p y/gd^X. Note that within a more general setup that includes lapse and shift functions, we would have to 
consider the question whether the locally lapsed time is replaced by a global "averaged time" that would involve an 
average over the lapse function. Here, the dust cosmology is already synchronous, so that this question does not arise. 
Note furthermore, that the building of averages is done in the inhomogeneous geometry. The averages functionally 
depend on the inhomogeneous metric, but this latter needs not to be specified. We may talk of a kinematical averaging 
that does not deform the geometry, i.e. that does not change the physical properties of the inhomogeneous space 
time. For other strategies, see [111, and references therein, as well as Section IV. 

• We generalize the kinematical laws of the standard model a) for the volume expansion (the Hamiltonian constraint in 
the ADM formulation of general relativity) and b) for the volume acceleration (Raychaudhuri's equation in the ADM 
formulation of general relativity) by dropping the symmetry assumption of local isotropy. The general equations are 
then volume-averaged, leading to the following general volume expansion and volume acceleration laws (for a volume 

1 /3 

scale factor, defined by (t) := (V^r>(i)/Vx)(ii)) j the overdot denotes partial time-derivative, which is the covariant 
time-derivative here) Q: 

3— --4^G(^^)I, + QD + A ; 3Hl + ^ ^S7:G{g)j,-lwv-lQv + I^, (2) 

where H-p denotes the domain dependent Hubble rate H-p = du/ax' = K is the trace of the extrinsic 

curvature of the embedding of the hypersurfaces into the space time, Kij, and A the cosmological constant. The 
kinematical backreaction Q-p is composed of averaged extrinsic curvature invariants, while Wp is an averaged intrinsic 
curvature invariant that describes the deviation of the average of the full (three-dimensional) Ricci scalar curvature 
TZ from a constant-curvature model, 

Qj,:=(k^-K)K^\ -l{K)l ; := (7^)^ - ^ . (3) 

\ / LJ O ^'P 

The kinematical backreaction Q-p can also be expressed in terms of kinematical invariants, where the extrinsic cur- 
vature is interpreted actively in terms of (minus) the expansion tensor: 

Qv:=l{{e')^-{Bfp)-2{a')^ , (4) 

where 9 is the local expansion rate and l/2<Tij(j'^^ is the squared rate of shear. Note that Hp is now defined 

as Hp — 1/3(0) J,. Q-p appears as a competition term between the averaged variance of the local expansion rates, 
(02)^ _ {0)1^ and the averaged square of the shear scalar (o'^)p on the domain under consideration. 

For a homogeneous domain the above backreaction terms Q-p and W-p, being covariantly defined and gauge in- 
variants in a perturbation theory on a homogeneous background solution, are zero. They encode the departure from 
homogeneity in a coordinate-independent way (ssl. l46|. 

The integrability conditions connecting the two Eqs. assuring that the expansion law is the integral of the 
acceleration law, read: 

(g) + 3Hp (g)^ = ; a:p''{alWp) + a^^a^^^Qp) = . (5) 



5 



While the mass conservation law for the dust is sufficient in the homogeneous case, there is a further equation 
connecting averaged intrinsic and extrinsic curvature invariants in the inhomogeneous case. The expressions in brackets 
are conformal invariants (for further details see [l8| . 

The interpretation of these average equations as generalized or evolving backgrounds [l8| , [4l| implies that the second 
conservation law describes an interaction between structure formation and background curvature. In the standard 
model this latter is absent and structures evolve independently of the background. This constant-curvature back- 
ground furnishes the only solution of ([5|) , in which structure formation decouples from the background (the expressions 
in brackets in the second conservation law are separately constant). Backreaction on such a fixed background decays 
in proportion to the square of the density and is unimportant in the Late Universe [ij, [3 ■ This degenerate case of 
a decoupled evolution explains the fact that in Newtonian and quasi-Newtonian models backreaction has no or little 
relevance [l^; in the Newtonian case Q, as well as in Newtonian [TJ, [38| and spatially flat, relativistic spherically 
symmetric dust solutions (s^ . Qp vanishes. In models with homogeneous geometry and with periodic boundary 
conditions imposed on the inhomogeneities on some scale, the backreaction term is globally zero and describes cosmic 
variance of the kinematical properties. 

Note here that, in general, a physical backgroimd "talks" with the fluctuations, and it is this coupling that gives 
rise to an instability of the constant-curvature backgrounds as we discuss below. The essential effect of backreaction 
models is not a large magnitude of Qx> , but a dynamical coupling of a nonvanishing Qx> to the averaged scalar curvature 
deviation Wv- This implies that the temporal behavior of the averaged curvature deviates from the behavior of a 
constant-curvature model. In concrete studies, as discussed further below, this turns out to be the major effect 
of backreaction, since it does not only change the kinematical properties of the cosmological model, but also the 
interpretation of observational data as we explain in Section IV. 



B. Scalar field emerging from geometrical inhomogeneities 

We rewrite the above set of spatially averaged equations together with their integrability conditions by appealing 

to the kinematical equations of the standard model, which will now be sourced by an effective perfect fluid energy- 
momentum tensor [10|: 

3!S=-47rG(pf^ + 3pSf) + A ; (6) 

3F|, - 2^ = STTGgff, + A ; (7) 

gf^ + 3Hv{g^^+pf^) - 0, (8) 

where the effective densities are defined as 

ef^ -.= {0)^ + 9^ ; := -_1^Q^ - -l^W^ ; 

p^^:^p., ; p., := -_1^Q^ + -l^W^ . (9) 

In this form the effective equations suggest themselves to interpret the extra fiuctuating sources in terms of a scalar 
field [ilH, which refer to the inhomogeneities in geometrical variables. Thus, we choose to consider the averaged 
model as a (scale-dependent) "standard model" with matter source evolving in a mean field of backreaction terms. 
This scalar field we call the morphon field, since it captures the morphological (integral-geometrical [Tsl ) signature 
of structure. (Note that in more general cases, involving lapse and shift functions, the structure of the scalar field 
theory suggested by the equations may no longer be a minimally coupled one.) We rewrite (l6| : 

qI^^\'^1 + Uv ; pI = ^\^1-Uv , (10) 

where e = +1 for a standard scalar field (with positive kinetic energy), and e = — 1 for a phantom scalar field (with 
negative kinetic energy; if e is negative, a "ghost" can formally arise on the level of an effective scalar field, although 
the underlying theory does not contain one; note also that there is no violation of energy conditions, since we have 
only dust matter). Thus, from the above equations, we obtain the following correspondence: 

^ -,Q-D = e^l-U-D ; -^yV-D = iUv ■ (11) 



SttG ' SttG 
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The correspondence ([TT]) recasts the integrability conditions ^ into a (scale-dependent) Klein-Gordon equation for 
$-p, and <i>-p 7^ 0: 

<i>v + 3Hv^v + e-^U{<i>v,{g)T,) = . (12) 

We appreciate that the deviation of the averaged scalar curvature from a constant-curvature model is directly pro- 
portional to the potential energy density of the scalar field. Averaged universe models obeying this set of equations 
follow, thus, a Friedniannian kinematics with a fundamental matter source, and an effective scalar field source that 
reflects the shape of spatial hypersurfaces and the shape of their embedding into spacetime. Given the potential in 
terms of variables of the averaged system, the evolution of these models is fixed (the governing equations are closed). 
This also potentially fixes coupling parameters, since all involved fields can be traced back to the initial value problem 
of general relativity. 

The morphon formulation of the backreaction problem opens a nice interpretation in terms of energies: a homo- 
geneous model, Qx> = (a necessary and sufficient condition to also drop the scale-dependence, if required on every 
scale), is characterized by the virial equilibrium condition: 

+ i^?ot= --^^ , Qv = ; i?S„ = £^?,Fp , E'^^, = -UvVv . (13) 

Deviations from homogeneity, Q-p ^ 0, thus invoke a non-equilibrium dynamics of the morphon in its potential that 
is dictated by the effective intrinsic curvature of the space in which the fluctuations evolve. Morphon energies are 
redistributed and can be assigned to the dark energies. Dependent on the signs of the backreaction terms (and a 
sign change may occur by looking at different spatial scales) the morphon can act as a scalar field model for Dark 
Matter, a quintessence model for Dark Energy, or it can even play the role of a classical inflaton, as we exemplify in 
the following subsection. (For the different interpretations of scalar fields see the review j25|, and for unified views 
the selection of papers [U, [H, [13, [S^l , and for scalar Dark Matter e.g. [3. [iol. [soj . 



C. Example: morphonic inflation 

Consider a tube of space time characterized by a gravitational field with no fundamental sources. The 4— Ricci 
curvature tensor vanishes everywhere, but not necessarily the 4— Weyl curvature tensor. Even if this classical vacuum 
space time is initially foliated into 3— Ricci flat hypersurfaces, this does not remain so in the dynamical evolution: 
such an initially prepared homogeneous state is unstable (a fact that we shall explain in the next section), and 
these hypersurfaces, if perturbed, will necessarily develop into inhomogeneous hypersurfaces featuring non-vanishing 
averaged curvature invariants, i.e. an intrinsic, on average negative curvature, and a compensating extrinsic curvature 
due to the embedding of the hypersurfaces into the Ricci-flat space time. Thus, in this picture, the space section will 
develop a morphonic scalar field that is driven by a Klein-Gordon dynamics and specified by the initial value of its 
self-interaction potential. While this instability is dynamical, the picture reminds us of the behavior of a fundamental 
inflaton, where the instability is created by an externally added potential. 

We specify initial data according to the analogy of the backreaction variables to the morphon field (Q^ = 
Qv{ti);Wiy = Wv{t^)): 

^---24b^- ' + ^ ' <f^-Mt.). (14) 

Interestingly, for a homogeneous initial state, Q-niti) — 0, the kinetic energy density of the morphon field is initially 
non-vanishing, and the Klein-Gordon dynamics drives the morphon into a stable fix point, an (assumed) existing 
minimum of the potential. However, the outcome does not depend much on the initial data for Qx>'. we could also 
start with inhomogeneous initial data, e.g. a cosmological constant that is mimicked by a particular morphon, in 
which case the initial kinetic energy density is zero. Curvature energy is thus converted into kinetic energy, driving 
the system into an accelerated expansion phase. The value of the potential is necessarily always positive, since the 
vacuum 3— space has negative intrinsic curvature. The infiationary mechanism is thus the same as the mechanism to 
create an accelerated expansion in a Quintessence phase. 

We realize the infiationary scenario [13 by closing the set of averaged equations with a potential of the generic 
Ginzburg-Landau form: 



(15) 
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FIG. 1: The Ginzburg-Landau potential in arbitrary units and the possible initial conditions as well as their physical meaning. 
All conditions possess some curvature yVx> < 0. The arrows schematically indicate the amplitude of the morphon's initial speed 

In the order of the points (from left to right): the first two points dominated by shear fluctuations (red, green) are obtained 
for Q^x> <0 <^ > 2{H^ + fc^,); the next points dominated by expansion fluctuations (blue, pink) for < 2{H^ + fc^), 
where the de Sitter-A equivalent case has a stiff morphon <I>^ = 0; the homogeneous case (last point, orange) is obtained for 

= 2{H\j^ + fc^). Figure from 



The position of the minimum $o and the amplitude Uq play different roles: the first one fixes the duration of inflation, 
while the second sets the size of the Hubble radius at which it happens. This potential has been extensively studied 
in the context of chaotic inflation (48| . The various initial conditions together with their interpretation in terms of 
geometrical properties of space are shown in Figure [TJ 

Combining this purely morphonic picture of inflation created from the Einstein vacuum with a fundamental scalar 
field, we can establish hybrid infiationary models with two scalar fields, one of them being the morphon that is always 
present in the case of inhomogeneous universe models. 



III. GLOBAL GRAVITATIONAL INSTABILITY OF THE STANDARD MODEL BACKGROUND 

A. The phase space of exact background states 

The space of possible states of an averaged cosmological model, or the space of "generalized backgrounds" has one 
dimension more than the space of possible homogeneous-isotropic solutions in the standard model framework. This 
can be seen by introducing adimensional "cosmological parameters". We divide the volume-averaged expansion law 
by the squared volume Hubble functional H-p '■= ap/ou introduced before. Then, the expansion law can be expressed 
as a sum of adimensional average characteristics: 



with: 



+ + nf: + + = 1 , (i6) 



O'D SttG / „\ _ . Q-D A . c>V "-^i 



Taking for simplicity the constant-curvature parameter and the curvature deviation into a single full curvature pa- 
rameter, + =: rt^, the generalized model offers a cosmic quartet of parameters. Furthermore, if we put A = 0, 
the expansion law defines, for each scale, a two-dimensional phase space of states. A one-dimensional subset of this 
phase space is formed by "backgrounds" with Friedmannian kinematics (see Figure [5]). 

We can analyze the fix points and their stability properties in the general dynamical system (l6| . (65| . The principal 
outcome of this study is that the standard zero-curvature model forms a saddle point, of particular interest are two 



8 




FIG. 2: Left: "Cosmic phase space" of the solutions of the averaged equations ("generaUzed backgrounds") in a plane s pan ned 
by the volume deceleration parameter := —d-D/ia-DH^) = 1/20^ + 2Qq — 57^ £^nd the matter density parameter [lq |. It 
represents a two-dimensional subspace { A = } of the full solution space that would include a cosmological constant. The 
segments are separated by particular exact scaling solutions of the full problem. We identify the following scaling solutions: 
all the scaling solutions are represented by straight lines passing through the Einstein-de Sitter model in the center of the 
diagram (1/2; 1). Models with "Friedmannian" kinematics, but with renormalized parameters form the line r = 1/3 (for details 
see [3, Appendix A). The line r = are models with no backreaction on which the parameter varies (scale-dependent 
"Friedmannian models"). Below the line r = in the "quintessence phase" we find effective models with subdominant shear 
fluctuations (Q-d positive, Qq negative). The line r = —1/3 mimics a "Friedmannian model" with cosmological constant. The 
line below r — —1/3 in the "phantom quintessence phase" represents the solution inferred from SNLS data (c/. [T^). and the 
point at {q^;Q,^) = ( — 1.03;0) locates the late-time attractor associated with this solution. Since we have no cosmological 
constant here, all expanding solutions in the subplane < drive the averaged variables away from the standard model 
featuring a backreact ion-driven volume acceleration of effectively isotropic cosmologies that are curvature-dominated at late 
times. 

Right: we show the evolution of phase space orbits (running away from the standard model) for a multiscale model that 
is explained further below; this model is partitioned into over-dense A4 (dotted) and under-dense £ (dashed) regions, their 
volume fraction being derived from A^'-body simulations (tH . They are shown here in the same plot for economic reasons and 
actually live in two different phase spaces corresponding to large scales for the f-regions and small scales for the A^-regions. 
denotes one of the regions M or £. 

instability sectors for the standard model, regarded as averaged state: firstly, perturbed homogeneous states are driven 
into a sector of highly isotropic, negative curvature and accelerated expanding "backgrounds" where backreaction thus 
mimics Dark Energy behavior over the domain T); secondly, perturbed homogeneous states are driven into a sector 
of highly anisotropic, positive curvature, collapsing and decelerated "backgrounds" where backreaction thus mimics 
Dark Matter behavior over the domain V. Concrete models show that the former happens on large scales, and the 
latter on the scales of galaxy surveys, and also on smaller scales. Thus, qualitatively, the instability sectors identified 
comply with the aim to trace the dark components back to physical properties, but they also agree with the expected 
properties of the structure: isotropic, accelerating states on large scales, and highly anisotropic structures on the 
filamentary distribution of superclusters. Moreover, the curvature properties also meet the expectations: on large 
scales the Universe is void-dominated and, hence, dominated by negative curvature, while on intermediate scales 
over-densities are abundant and are characterized by positive curvature. 



B. Dark Energy and Dark Matter hidden in the geometry of space 

The fact that the standard model is globally unstable in the phase space of averaged states, and the fact that 
the instability sectors lie in the right corners to explain Dark Energy and Dark Matter behavior, are botlr strong 
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qualitative arguments to expect that the conservative explanation of the dark energies through morphon energies is 
valuable. The underlying mechanism is indeed based on the fundamental existence of the relation between geometrical 
curvatures and sources dictated by Einstein's equations. 

Whether this mechanism is sufficient in a quantitative sense is to date still an open issue. The difficulty to 
construct quantitative models is to be seen in the need for non-standard tools, for example perturbation theory on 
a fixed reference background should be replaced by a fluctuation theory on evolving backgrounds that captures the 
average over the fluctuations. The question whether perturbations are small can only be answered if we know with 
respect to which background they are small. Furthermore, since backreaction affects the geometry, it will change the 
interpretation of observational data, a problem that is intimately related to the generalization of the cosmological 
model, and to which we shall come below. 

Before, we shall in the next section explain the identified mechanism by discussing some physical properties of 
structure formation and its relation to the interpretation of geometrical curvature invariants and how they mimic the 
dark sources. We here touch on a deeper problem: backreaction effects account for both. Dark Energy and Dark 
Matter, simultaneously. Whether, on a given domain, or on an ensemble of domains on a given scale, the morphon 
mimics Dark Energy or Dark Matter behavior, changes as a function of time and as a function of scale. Moreover, 
the small-scale contribution to e.g. a Dark Matter behavior requires more sophisticated relativistic models than the 
dust model used throughout here (e.g. (H^l). Considering rotation curves of galaxy halos, missing gravitational 
sources in clusters or missing sources on cosmological scales always needs different modeling strategies. We try in the 
following to provide a first step of disentangeling Dark Energy and Dark Matter behavior by explicitly constructing 
an effective multiscale cosmological model. 

C. Multiscale cosmology and structure— emerging volume acceleration 

Contrary to the standard model, where a homogeneous background is used as a standard of reference for the 
expansion history of the Universe, a background constructed as the average over fluctuating fields introduces a subtle 
element: while a homogeneous geometry can be characterized locally, an average is nonlocal, since it is determined 
by the inhomogeneities inside, but also outside the averaging domain, reflecting the nonlocal nature of gravitation. 
Furthermore, an average incorporates correlations of the local fields. It is this latter which is the key-driver of a 
repulsiv "effective pressure" that arises in the averaged models. 

This "effective pressure" provides the reason why backreaction can produce a volume-accelerating component despite 
the decelerating nature of the general local acceleration law. This can be seen easily by comparing the local and the 
volume-averaged Raychaudhuri equation (for vanishing vorticity and pressure that both would also act accelerating 
on the local level, but only on small scales): 

^ = A-47rG£< + 2II-l2 ; (61) = A- 47rG(^?)p + 2(11) J,- (I)^ , (18) 

where we defined the principal scalar invariants of the expansion tensor 9^^ , 211 := 2/36^ — 2tT^ and I := 9. 

Clearly, by shrinking the domain to a point, both equations agree. However, evaluating the local and averaged 
invariants, 

2II-l2 = -i02„2ff2 . 

2 (11)^ - (I)^ = ! ((^ - Wv)')^ - 2 ((a - (a)^)2)^ - i {9)1 2 {a)l , (19) 

gives rise to two additional, positive-definite fluctuation terms, where that for the averaged expansion variance enters 
with a positive sign. Thus, the time-derivative of a (on some spatial domain 2?) averaged expansion may be positive 
despite the fact that the time-derivative of the expansion at all points in V is negative. 

In concrete models this variance is the source of a possible large-scale volume-acceleration that would be assigned 
to Dark Energy in the standard model, while the averaged shear fluctuations mimic an attractive source that would be 
missing as Dark Matter in the standard model on cosmological scales. Both terms are competing in the kinematical 
backreaction Qx>. Since backreaction depends on scale, it may act in both ways. 

We can go one step further and make the scale-dependence explicit by introducing a union of disjoint over-dense 
regions M and a union of disjoint under-dense regions £, which both make up the total (honiogeneity-scale) region 
v. The averaged equations can be split accordingly yielding for the kinematical backreaction [131 , |7l| : 



(20) 
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FIG. 3: Plot of the evolution of Q-d and {TZ)jj in terms of the global scale factor o-d. For comparison a line with a simple 
^ scaling is added, that corresponds to the leading mode in second-order perturbation theory, and one that is constant. 
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Qjy and (7?.)^, are normalized by — 6i?|). so that the values at the initial time represent fJg' and 0,^ . We appreciate that the 
backreaction terms feature an approximate cosmological constant behavior on the homogeneity scale despite the assumption of 
a^^-scaling on the partitioned domains T — M,£. Physically, this result can be attributed to the expansion variance between 
the subdomains and, hence, this latter is identified as the key effect to produce a global Dark Energy-like behavior of the 
backreaction terms. Figure from [Tlj . 



where \m '■— / denotes the volume-fraction of the over-dense regions compared to the volume of the region 
T). In a Gaussian random field this fraction would be 0, 5 and would gradually drop in a typical structure formation 
scenario that clumps matter into small volumes and that features voids that gradually dominate the volume in the 
course of structure formation. 

Ignoring for simplicity the individual backreaction terms on the partitioned domains, the total backreaction features 
a positive-definite term that describes the variance between the different expansion histories of over- and under-dense 
regions. It is this term that generates a Dark Energy behavior over the domain V (see also [H^ for a model by Rasanen, 
and (Zsl . Ft^ - It^ for Wiltshire's model that is based on this term only, but includes a phenomenlogical lapse function 
to account for different histories in Ai and £ regions that, this latter, we cannot implement in the synchronous 
foliation of a multiscale dust model). If we model non-zero individual backreaction terms by an extrapolation of the 
leading perturbative mode in second-order perturbation theory (46l . |47| that also corresponds to the leading order 
in a Newtonian non-perturbative model |11| . then we even produce a cosmological constant behavior over I?, see 
Figure [21 In other words, the fact that, physically, over-dense regions tend to be gravitationally bound, i.e. do not 
partake significantly in the global expansion, already produces a large-scale "kinematical pressure" as a source of 
volume acceleration. A homogeneous background simply cannot account for this difference. 
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IV. INHOMOGENEOUS AND EFFECTIVE METRICS 
A. Some notes on relativistic perturbation theories 

Consider a spatial metric form g with coefficients gij in an exact (co-tangential) basis dX^ (E}dX^ . We can write any 
metric as a quadratic form of deformation one-forms, g = 5ab'n°'®'ft ^ i-S- in terms of coefficients, = SabV'^iV^j- Now, 
such a metric form is homogeneous, i.e. its Ricci tensor vanishes everywhere, if there exist functions /°, such that the 
one-forms can be written as exact forms, = d/". In other words, if we can find a coordinate transformation = 
f°-^^{X^,t) that transforms the Euclidean metric coefficients in a new basis, dx'^(E)dx^ , Sijdx^dx^ = Sabf^^if^^jdX^dX^, 
with a vertical slash denoting partial spatial derivative, into the metric coefficients gij, then these latter are just a 
rewriting of the homogeneous space. Given this remark, any perturbation theory that features metric forms of the 
integrable form, does not describe relativistic inhomogeneities; metric coefficients of the form gij = Sabf^^if^^j describe 
Newtonian (Lagrangian) perturbations on a flat background space. A truely relativistic perturbation theory deforms 
the background geometry, in other words, the perturbations live in a perturbed space, not on a reference background. 
This remark also shows that relativistic perturbation terms can never contain full divergences, since this latter needs 
integrable one-form fields. 

In light of these introductory remarks, an inhomogeneous relativistic metric produces curvature that, if volume- 
averaged on some domain, does not obey a conservation law (as can be explicitly seen in the coupling equation to 
the fluctuations ([5])) in the sense that it would always average out to zero; for details on curvature estimates see 
[TtI). This fact in itself shows the existence of a dynamical evolution of an averaged curvature, as soon as structures 
form. On the contrary, standard perturbation theory formulated on a fixed flat background is constructed such that 
the averages always vanish on the background, demonstrating the limited nature of results obtained by standard 
perturbative arguments. 

Another perturbative argument aims to justify the validity of the homogeneous geometry, even down to the scales of 
neutron stars (37| . since perturbations of the metric remain small with respect to the fiat background. This argument 
does not contradict the existence of a large backreaction effect, since these latter depend on first and second derivatives 
of the metric [l^, [l^l, Also, the perturbations are considered on a flat background that does not interact with 
structure. As we explained in detail, the perturbations may be small on a different (physical) background, in which 
case a perturbation may already live in a background with strong curvature (a zero-order effect). It is therefore 
not fruitful to argue against the relevance of backreaction within standard limited schemes, but rather an effort to 
generalize perturbation theory is needed. 

B. Template metrics and effective distances 

For the construction of an effective cosmological evolution model, as outlined above, a metric needs not be specified. 
The need for the construction of an effective metric in these models arises, since measured redshifts have to be 
interpreted in terms of distances along the light cone. Given an explicit, generic and realistic, inhomogeneous metric, 
the need for the construction of effective metrics does not arise. Also, if we succeed to understand the evolution of 
light cone averages in relation to distances, then also here an explicit metric will not be needed (goI. |6T|. (2l| . 

The idea of an effective cosmological metric comes from the "fitting problem", that has been particularly emphasized 
by George Ellis already in the 70 's [2^. The observation was that an inhomogeneous metric does not average out to 
a homogeneous metric that forms a solution of general relativity. Not only the nonlinearity of the theory, but also 
simple arguments of a non-comniutativity (28| between evolution equations and the averaging operation, give rise to 
the need to find a "best-fit", we may call it "template" geometry, that inherits homogeneity and (almost-)isotropy 
on the large scales and, at the same time, incorporates the inhomogeneous structure "on average" (see also the early 
practical implementations of this problem [3l|,[33, [13, (3^). 

For the solution of the fitting problem various strategies have been proposed (see and references therein). One 
strategy, that allows to explicitly perform a "smoothing" of an inhomogeneous metric into a constant-curvature metric 
at one instant of time, is based on Ricci-flow theory: one notices that a smoothing operation of metrical properties 
can be put into practice by a rescaling of the metric in the direction of its Ricci curvature. The scaling equations 
for realizing this are well-studied, and the rescaling flow results in a constant-curvature metric that carries "dressed" 
cosmological variables [l^l , [I^l ■ These incorporate intrinsic curvature backreaction terms describing the difference to 
the "bare" cosmological parameters as they are obtained through kinematical averaging. 
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FIG. 4: Left: Comparison of the luminosity distances of the multiscale models investigated in |71| and based on the template 
metric of [43 |. with the one of a flat ACDM model with h — 0.7 and Qrn ~ 0.27. On top a model where we force the volume 
scale factor a-D to follow the ACDM evolution. Despite this assumption, the changing curvature affects the luminosity distance. 
The luminosity distances in these models show a significant feature at a redshift of around 1, when compared with the best fit 
ACDM model, which may be looked for in the SN data. The curve below is a model with a^^-scaling. For comparison we also 
included the luminosity distance of the best fit model of [43 |. Because of a different Hubble rate of /i = 0.7854 it lies below 
the others from the beginning. This model does not significantly show the distinct feature of the other two models around a 
redshift of 1, due to the assumption of a sing le-scale cosmology. 

Right: Values of Clarkson's C-function [23l | for the best fit model of [i^ l (top), the model where the scale-factor is forced 
to follow the ACDM evolution (middle), and the model with a^^-scaling (bottom). Recall that, for every Friedmann model, 
C{z) vanishes exactly on all scales and for all redshifts. For the inhomogeneous models shown in the plot, this function has 
a minimum which may serve as observational evidence for the effective cosmologies, as proposed in [4a |. As both multiscale 
models shoWj_2t is not even necessary to measure derivatives of distance, since the feature is already present in the distance. 
Figure from [71| . 



C. Reinterpretation of observational data 

The standard method of interpreting observations is to construct the light cone ds^ — from the line-element 
ds'^ = -dt'^ + gf°'^dX^dX^ , where the coefficients gf°"^ are given in the form of a constant-curvature (FLRW) metric, 
and then to calculate the luminosity distance ^^(z) in this metric for a given observed redshift z. Assuming this metric 
for the inhomogeneous Universe implies the conjecture that the FLRW metric is the correct "template" of an effective 
cosmological metric. However, the integrated exact equations (the integral properties of a general inhomogeneous 
model) are not compatible with this metric, simply because the averaged curvature is assumed to be of the form 
(TV) J, = Qka"^ on all scales. Improving the metric template slightly, by replacing the global scale factor a{t) through 
the volume-scale factor ax>{t) and the integration constant k through the domain-dependent integration constant 
kx), renders this metric implicitly scale-dependent (53| . As we explained, this is not enough since the averaged 
curvature couples to the inhomogeneities and in general deviates from the a^^-behavior. What we can do as a first 
approximation, and this would render the metric compatible with the kinematical average properties, is to introduce 
the exact averaged curvature in place of the constant curvature in this metric form (43| . 

The resulting effective space time metric consists of a synchronous foliation of constant-curvature metrics that 
are, however, parametrized by the exact integral properties of the inhomogeneous curvature, thus they "repair" the 
standard template metric as for the evolution properties of spatial variables. Such a construction can be motivated 
by Ricci-flow smoothing, that guarantees the existence of smoothed-out constant curvature sections at one instant of 
time, and by assuming that the intrinsic backreaction terms are subdominant, so that we can parametrize the metric 
by "bare" kinematical averages. To stack these hypersurfaces together introduces, however, an inhomogeneous light 
cone structure [5l|, [H^. Ideally, one would wish to smooth the light cone too, which is also possible by employing 
Ricci flow techniques. Improving this first approach to a template metric is needed and this is work in progress. 

The result of employing the improved template metric described above is a change in the luminosity distance 
that would alter the interpretation of all observational data formerly based on FLRW distances. Examples for the 
multiscale models investigated in {Tlj are presented in Figure^ Although this investigation certainly needs refinement, 
we already appreciate a signature of the different curvature evolution that furnishes a clearcut prediction for future 
observations (see [i^ for details). 
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